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The present status of calculations of four fermion final states is reviewed. Higher order problems arising there 
are pointed out. Special attention is paid to results obtained in the semi-analytical approach and their limitations 
and perspectives. 



> 

(N 

\o 

(N 

\o 

O 

On 

Oh! 
CD 



X 



1. Introduction 

Four fermion final states are already observed 
at LEP1 [ 0. They are the main final states to 
be investigated at LEP 2 [ ^] . The two LEP runs 
planned this year 

June/ July at ^/s = 161GeV 
September/October at ^/s — \7AGeV 

will give the first opportunity to compare theoret- 
ical predictions with the experimental data at the 
highest energies ever reached in e + e~ collisions. 

The production of two jermions in the final 
state is quite simple at the Born level, ft pro- 
ceeds through annihilation of gauge bosons in the 
s channel, figure la. For Bhabha scattering, the 
gauge boson exchange is also possible in the t 
channel, figure lb. 

Four fermion final states can be created only 
in higher order processes. They are described by 
sets of Feynman diagrams, which have a much 
richer topology. Two diagrams, figure lc and Id 
are essentially new containing the non-abelian in- 
teraction of gauge bosons. The diagrams le and 
If can be interpreted as the radiation of a fermion 
pair from the initial and final state of the dia- 
gram la. Similarly figures lg and lh show pair 
radiation from the diagram lb. The correspond- 
ing diagrams with Higgs exchange are not shown 
in figure 1 . They can be obtained by replacements 
of a gauge boson with the Higgs at the appropri- 
ate places. 

The cross sections of four fermion final states 
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Annihilation 
Fig. la 



Exchange 
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Annihilation 
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Annihilation Conversion 
Fig. le Fig. If 



Bremsstrahlung 
Fig. lg 



Bremsstrahlung 
Fig. lh 



Multiperipheral 
Fig. li 



Fig. 1: Topologies of the lowest order Feynman di- 
agrams for e + e~ — » // and e + e~ — > /1/2/3/4- The 
dashed lines represent any allowed gauge boson. 



can be of the order of those of fermion pair pro- 
duction due to large factors, which can arise and 
compensate the additional a 2 . 

The diagrams in figures lc and If contain two 
resonating gauge bosons, which tend to be on 
mass-shell. They are called signal diagrams for 
W and Z pair production. Any resonating gauge 
boson V enhances the cross section by roughly a 
factor 



37 for V = W or Z. 
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If the resonating boson is a photon going to 
//, the enhancement factor (without cuts) is 
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ln[s/(m(/) + m(f)) 2 ]. If the photon couples 
to light fermions, we recover the infrared diver- 
gence due to soft photons. We assume that both 
fermions /, / are seen in the detector. This im- 
plies that the virtual photon cannot be arbitrarly 
soft. The diagrams shown in figures Id, le, lg 
and lh are enhanced by only one resonating gauge 
boson. 

Diagram le is unique because it is the only 
one, which allows gluon exchange at the tree level. 
This gives an enhancement by a factor 



14 2 to 25 2 



(2) 



depending on the scale q 2 of a s . 

All diagrams with gauge boson exchange in the 
t channel (figures lb, Id, lg, lh and li) contain 
collincar singularities. Although they are regular- 
ized by the non-zero electron mass, an enhance- 
ment factor 
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remains for every collinear photon in the total 
cross section. This large factor can be reduced by 
a cut on the angle between the outgoing e + (e~) 
and the beam axis, 

Let us be more explicit and estimate the total 
cross sections for e + e _ — > W + W~ and e + e~ — > 
e + e~ ff compared to e + e~ — > // by collecting 
only large enhancement factors. 

W pair production has a double resonant 
enhancement. A further enhancement comes 
from the counting of the possible final states. 
Fermion pair production results in 5 different 
quark flavours with three different colors each or 
in 3 different lepton or neutrino pairs. This adds 
up to 21 different final states. A W can decay in 
2 different quark pairs of 3 different colors each 
or in 3 different lepton pairs. This gives 9 differ- 
ent final states. We have yet to take into account 
the combinatorical factor v2j which enhances the 
W /1/2 couplings relative to the Zff couplings. 
This is also the reason, why T 'z/^w ~ 21/(2-9) 
and not 21/9. Combined, we get the following ad- 
ditional factor for W pair production compared 
to fermion pair production: 
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After dividing out the threshold factor j3 = 
yl — 4Af^/s for W pair production, this crude 
estimate is correct within a factor 1.5 for 2Mw < 
\/s < 500 GeV. For very high energies, the cross 
section for W pair production becomes larger 
than the estimate (Q) due to terms proportional to 
ln(s/M|p). The production of two intermediate Z 
bosons is suppressed by the effective symmetriza- 
tion arising compared to W pair production. 

The cross section of the reaction e + e~ — * 
e + e~// is dominated by diagrams shown in fig- 
ure li. All three propagators in the t channel can 
become singular simultaneously leading to huge 
enhancement factors. These factors can be esti- 
mated starting from the description of the pro- 
cess as a scattering of photons with Weizsacker- 
Williams spectrum [ ||) . Keeping only the leading 
terms, one gets 



er(e + e -> e+e //) 
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For / = b, we get at y/s = 200 GeV 

a(e + e~ — > e + e~bb) 
w a QED (e + e- bb) ■ 3 ^ 2.5pb. (6) 

This is less than a factor 2 off from the exact 
result, see figure 13 in [ |J. 

All four fermion final states can be classified 
according to the topologies of Feynman diagrams 
entering the process. The number of the Feyn- 
man diagrams involved in different processes with 
charged current exchange is shown in table 1, 
while table 2 contains final states, which are pro- 
duced by the exchange of neutral gauge bosons 
only. Potential diagrams with Higgs exchange 
are not counted. Diagrams coming from different 
contributions of quark mixing are counted. Final 
states indicated in bold contain only the topolo- 
gies of figures le and If. In the case of charged 
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currents, the topology lc has to be added. These 
charged current processes are usually referred as 
CC9, CC10 or CC11 processes. Numbers in ital- 
ics indicate final states, which belong to table 1 
and table 2. The CC43 and CC19 (identical to 
NC43 and NC19) processes are described only by 
the topologies of the bold processes. However, 
here new structures in the interferences appear. 
The CC56 (NC56) process involve all topologies. 
All remaining processes contain electrons or elec- 
tron neutrinos in the final state allowing gauge 
boson exchange in the t channel. The correspond- 
ing numbers are written in roman. They contain 
all topologies simultaneously, except diagram Id, 
which demands at least one electron neutrino in 
the final state. Final states with identical final 
fermions are given in typewriter. They require 
additional contributions to insure Pauli antisym- 
metrization. 
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Tab. 1: Number of Feynman diagrams for CC type 
final states, see [\2a^. 
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Tab. 2: Number of Feynman diagrams for NC type 
final states, see [\24y- 

2. Higher order corrections 

Four fermion final states can be produced with 
cross sections of several pb. With an integrated 
luminosity of 500pfe _1 per year, the expected ex- 
perimental precision for some cross sections is 



about 1%. It has to be met by theoretical predic- 
tions, which have errors of 0.5% or less. Hence, 
higher order corrections have to be included. 

2.1. QED corrections 

The complete QED corrections to four fermion 
final states are not known. One could simplify 
the problem for W pair production considering 
corrections to the signal diagrams figures lc and 
If only. As an important step, this was done in 
the case of on-shell W's [ §. For off-shell W% 
the complete correction is not known. 

In the following, the Coulomb singularity and 
initial state corrections are mentioned. 

Coulomb singularity 

The Coulomb singularity [ |E| arises from long 
range electromagnetic interactions between the 
produced massive charged particles. For W pair 
production, we get the correction 



Coul n 
a = a 
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to the Born cross section cr°, which diverges near 
threshold where the velocity of the VF's j3 = 
yl — 4M^ / s approaches zero. It indicates that 
perturbation theory is not applicable in this re- 
gion. Fortunately, the non-zero width T\y and a 
slight off-shell production of the VT's regularize 
the singularities in equation (Q). The numerical 
effect can exceed 6% of the total cross section 
near threshold [ §|, 0]- Therefore, the Coulomb 
correction has to be inlcuded. 

Formula (0) also applies for QED and QCD cor- 
rections to pair production of massive fermions as 
it appears in the diagrams ld-lh. This is well 
known from final state corrections to the dia- 
gram la at LEP 1 [ |8| . The problems can be cured 
by a calculation in the limit of massless fermions 
or by a cut on the invariant mass of the heavy 
fermion pairs. Such a cut is desirable for quark 
pairs in any case to avoid non-perturbative bound 
state regions. 

Initial state corrections 

Initial state QED corrections to off-shell W pair 
production are calculated in [ |^] . They reach sev- 
eral % near the WW threshold. In such calcu- 
lations a problem arises because it is not clear 
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how to separate weak corrections from QED cor- 
rections and initial state QED corrections in a 
unique and gauge invariant way. This separa- 
tion problem is solved by the current splitting 
technique [ 0], in which the chargeless neutrino 
exchanged in the t channel is divided into two 
charge flows of opposite sign. Now the charge 
flows of the initial and final states are separated, 
and gauge invariance is ensured as it is in the case 
of Z pair production. The resulting initial state 
QED corrections can be split into universal con- 
tributions, which are described by the same flux 
function as known from annihilation diagrams in 
the case of LEP 1 physics, and into non-universal 
contributions depending on the particular pro- 
cess. The non-universal contributions to off-shell 
W and Z production are numerically small for 
LEP 2 energies. They are suppressed by a factor 
S1S2/ ' s 2 , where si^ are the invariant energy flows 
through the JU's (Z's) [ |, |§. The smallness of 
the non-universal corrections cannot be taken for 
granted, but has to be proven for any particular 
process. 

The universal corrections lead to handy formu- 
lae for initial state corrections to cross sections of 
W and Z pair production, 

a is R{s) = r u^ j j 

Js>- S 

Alternatively, initial state radiation can be taken 
into account by the structure function approach [ 
[TlJl . It assumes that both colliding photons have 
energies degraded by radiated collinear photons, 

o* s «{8) ^ (9) 

= / dx% / dx2-D(xi, s)D(x2, s)o~°(xiX2s). 

The flux function (FF) p(s'/s) and the structure 
function (SF) D{x, s) contain information about 
real and virtual corrections. Further details can 
be found, for example, in [Q] or [ [l2|| . 

2.2. Electroweak corrections 

Complete one-loop corrections to four fermion 
processes would be the best task to get more ac- 
curate descriptions of these cross sections. Un- 
fortunately, this calculation is very complex [ |l3| ] 
and done only for on-shell Ws [ ||. 



In this paragraph, only electroweak corrections 
connected with the finite widths of gauge bosons 
are mentioned. 

At lowest order, a particle V has no width, 
i.e. its propagator Dy is proportional to Dy ~ 
[q 2 — My]" 1 , where q 2 is the momentum transfer 
through the particle of mass My. 

Massive gauge bosons are unstable having a 
non-negligible decay width. The required preci- 
sion at the energy range of LEP 2 demands the 
inclusion of a finite width. Although the inclu- 
sion of a width for particles in the t channel is 
not physical, the simplest possibility is the in- 
clusion of the constant on-shell decay width Ty 
everywhere in the propagator, Dy ~ [q 2 — M v + 
iMyTy}- 1 . 

The next step of accuracy, which cures this de- 
ficiency, is to take into account the natural energy 
dependence of the width arising from self-energy 
insertions into the boson propagator. Present and 
future experiments are sensitive to the difference 
of these handlings of the width [ [l4|, [l5| , even at 
the limit of very high energies. 

A finite width in the propagator implies the 
partial inclusion of contributions, which are of 
higher order in perturbation theory. Unfortu- 
nately, the symmetries of a theory (as gauge in- 
variance) are only respected in a fixed order of 
perturbation theory. Hence, the inclusion of a 
non-zero width leads to gauge dependent cross 
sections. Although the gauge dependent terms 
are of higher order in perturbation theory, they 
can be enhanced by large kinematical factors as 
s/m 2 in the reaction e + e~ — * e~D e ud, see [ 
Therefore, finite widths must be included with 
care. Different approaches for the restoration 
of gauge invariance have been proposed, among 
them the Laurent expansion of the matrix ele- 
ment and the inclusion of projection operators [ 
[ijj or the inclusion of certain higher order contri- 
butions [ m. 

2.3. QCD corrections 

QCD corrections give sizeable contributions to 
distributions and cross sections. They are known 
from LEP 1 for the production of pairs of (heavy) 
quarks [ || . They are under investigation for four 
quarks in the final state [ |l^] . 
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QCD corrections enter the width of the Z, W 
and the Higgs. For the case of rV, QCD correc- 
tions can reach several %, see [ |2^]. Most of the 
QCD corrections to the Higgs width can be ab- 
sorbed into the running quark masses evaluated 
at the scale of the Higgs mass. They lead to values 
for m b {M H ) = 2.9 GeV and m c (M H ) = 0.6 GeV 
[ |l| , which deviate substantially from the corre- 
sponding pole masses. 

2.4. Theoretical uncertainties 

The theoretical uncertainties are extensively 
discussed in the different contributions to [ ||. 
Here, we only give some of them in a telegraphic 
style. 

QED corrections: 

• The initial state QED corrections to all Born 
cross sections presented in [ || are calculated by 
the FF or by the SF approach. The validity of 
these approximations is proven only for W and 
Z pair production and for annihilation diagrams. 
For other processes it is not known, whether the 
non-universal corrections are small enough to be 
neglected. Even for W and Z pair production, 
they reach 0.3% at LEP2 energies and rise above 
1% at ITeV 

• The FF is derived after an average over all phase 
space variables except the two invariant energy 
flows s 12 through the W's or Z's. For distribu- 
tions or cuts in parameters different from s\ and 
S2 (i.e. angular cuts, cuts on energies of single 
particles) , it depends in general on the additional 
parameters. This is known from LEP 1, where the 
angular-dependent FF for annihilation diagrams 
is derived [ ^2|. The dependence on the distri- 
bution parameter disappears only for soft photon 
radiation. Therefore, the theoretical errors of ob- 
servables, which favor hard photon radiation can 
become very large. Fortunately, these observables 
have small cross sections because QED favors soft 
photon radiation. 

• We now assume that the FF and SF approaches 
work perfectly to all orders of perturbation the- 
ory. We can then use our knowledge about these 
functions beyond the leading order and study the 
changes of cross sections. This relative change 
reaches 0.5% [ |2(| giving an idea about the effect 



of higher order QED corrections. 
QCD corrections: 

• A calculation of complete QCD corrections to 
four fermion processes is not done. The result- 
ing uncertainty can be estimated in the follow- 
ing simple example: W pair production with 
four quarks in the final state involves diagram le 
with gluon exchange at the Born level. The di- 
agram favors soft gluons. The minimal momen- 
tum transfer through the gluon is usually defined 
by cuts. The scale q of a s (q 2 ) in these contri- 
butions is not known exactly. An uncertainty 
of 20% in the scale q leads to an uncertainty of 
10% in a 2 (q 2 ) for q w 10 GeV. If the diagram le 
would contribute 10% to the total cross section, 
this scale uncertainty would transform to an er- 
ror of 1% in the total cross section. If the cuts 
allow smaller q 2 , a s (q 2 ) and its uncertainty be- 
come larger. This simple example illustrates the 
dependence of QCD corrections on kinematical 
cuts. 

• The running quark masses in Higgs decay de- 
pend on a s (Mz). The resulting uncertainties are 
shown in table 3 in [ [2l[]. They reach several %. 
They transform into uncertainties of predictions 
for production cross sections of Higgs Bosons be- 
ing proportional to the square of these masses. 

• The interface between partonic cross sections 
and hadronization procedures and the hadroniza- 
tion procedures itself introduce theoretical uncer- 
tainties. 

Weak corrections 

Gauge invariance can be restored in calculations 
with finite widths by different methods. The aris- 
ing numerical difference is studied for the process 
e~e + — > e~P e ud in the second reference of [ 18 



and found to be about 0.5%. The disagreement 
between the different methods is expected to be 
larger near threshold [ [l3| . 

Final remark 

The estimation of theoretical errors is a highly 
subjective task. These errors can be reduced only 
by future higher order calculations. Remember- 
ing that not all theoretical errors are equally im- 
portant in all observables, we conclude that the 
present theoretical accuracy probably meets the 
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precision expected in the first year at LEP2. It 
has certainly to be improved in the future. 

3. Codes and algorithms 

3.1. Existing codes 

Event generators for four fermion final states 
are described in detail in volume 2 of reference [ 
[|. They use different algorithms for the evalua- 
tion of the squared matrix elements, phase space 
integration and mapping of singularities. Numer- 
ical comparisons between the programs gave a 
nice agreement. 

It follows a list of the available codes and its 
authors. 

Monte Carlo Programs: 

• ALPHA (F. Caravaglios, M. Moretti) 

• CompHep 3.0 (E. Boos, et al.) 

• ERATO (C.G. Papadopoulos) 

• EXCALIBUR (F.A. Berends, R. Kleiss, 
R. Pittau) 

• grc4f 1.0 (J. Fujimoto, et al.) 

• KORALW 1.03 (M. Skrzypek, S. Jadach, 
W. Placzek, Z. Was) 

• LEPWW (F.C. Erne) 

• LPWW02 (R. Miquel, M. Schmitt) 

• PYTHIA 5.719 / JETSET 7.4 (T. Sjostrand) 

• WOPPER1.4 (H. Anlauf, T. Ohl) 

• WPHACT (E. Accomando, A. Ballestrero) 

• WWF 2.2 (G.F. van Oldenborgh) 

• WWGENPV/HIGGSPV (G. Montagna, 
O. Nicrosini, F. Piccinini) 

Other codes: 

• GENTLE/4fan (D. Bardin, D. Lchner, A. Leike, 
T. Riemann), Semi-analytical code 

• WTO (G. Passarino), multi-dimensional 
deterministic integration 

3.2. Comparison between the Monte Carlo 
the and semi-analytical approach 

Many different topologies of Feynman diagrams 
are involved in calculations of four fermion pro- 
cesses leading to a rather singular matrix element. 
It has to be integrated over the 8-dimensional 
phase space. In the case of absence of transver- 
sal beam polarization, one degree of freedom, the 
rotation around the beam axis, is trivial. We 



are left with a 7-dimensional phase space inte- 
gration. In the semi-analytical (SA) approach, 
most of the phase space integrations are done an- 
alytically. The remaining integrals are calculated 
numerically. In the Monte Carlo (MC) approach, 
all integrals are taken numerically. 

As shown in section 2, kinematical cuts are 
needed to reject threshold regions of massive 
fermion pairs because they are not described by 
perturbation theory. Additional angular and en- 
ergy cuts are required to ensure the detection of 
four final particles in a real detector. Finally, a 
real detector has wholes, which are planned (i.e. 
for cables and support) and those which arise 
spontaneously (i.e. damaged and dead sectors). 
These individual properties have to be inlcuded 
into a data analysis. There is an additional phys- 
ical motivation for cuts to separate interesting 
events from the background. 

Having in mind these requirements, it is obvi- 
ous that one needs a MC event generator. It is 
most flexible in kinematical cuts and can be con- 
nected with detector simulations and hadroniza- 
tion procedures. 

On the other hand, fits require fast and ac- 
curate programs. The overall accuracy (which 
consists of theoretical uncertainties and numer- 
ical uncertainties) must not be larger than half of 
the experimental error. Remembering the present 
theoretical errors, it follows that the numerical 
uncertainty should be considerably smaller than 
the experimental error. Demanding that a fit to 
future data should be finished rather in 1 to 10 
days and not in 10 to 100 days, and assuming 100 
to 1000 cross section calculations for one fit, one 
can estimate the calculation speed and accuracy 
of a code, which is useful for fits: Such a code 
should calculate one cross section with an numer- 
ical accuracy of < 0.1% in less than 15 minutes 
on computers already available today. The higher 
speed of future computers will probably be eaten 
by the implementation of future theoretical re- 
sults needed to decrease the theoretical errors. 

These requirements can be met by a SA code. 
In the simplest case of the NC32 processes (all 
bold final states in table 2), the typical calcu- 
lation time is 30s with a numerical accuracy of 
0.1% and with QED corrections included by the 
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FF approach. This time can vary by a factor 3 
depending on the cuts and on the final state. The 
calculation time for CC11 processes (all bold fi- 
nal states in table 1) is about 3 times longer due 
to the longer analytical formulae involved. The 
codes of GENTLE/4fan distributed during the 
LEP 2 workshop are much slower. This will be 
improved with an update at the end of this year. 

SA codes are not as flexible as MC's allowing 
only the calculations of a limited number of distri- 
butions. However, they establish a benchmark for 
MC's because SA codes have an inherently much 
smaller numerical error. Unfortunately, the SA 
method can fail for certain processes. For details, 
see the next section. 

The calculation time of SA and MC programs 
depends on the numerical accuracy e required. In 
the SA approach, we have an additional depen- 
dence on the dimension d of the numerical inte- 
gration and on the order r of the numerical inte- 
gration algorithm: 



MC: t 
SA: t~ 



e 

c -d/r 



(10) 

A few remarks are in order inspecting relations 
©: 

• In practice, the calculation time for MC's de- 
pends on d because the integrand cannot be 
mapped to a completely flat function. 

• High accuracy demands large calculation times 
in the MC approach. The difference between 
e = 1% and e = 0.1% has to be paid by a factor 
100 in computer time. 

• In the SA approach, we have r = 5 for Simp- 
son's rule. More sophisticated algorithms with 
larger r don't necessarily lead to shorter calcula- 
tion times because their errors are proportional 
to higher derivatives of the integrand. 

• SA codes can easily achieve a high accuracy 
They are very fast for small d. They fail for large 
d. In practice, the calculation time increases by 
more than a factor 10 for every additional integra- 
tion. We have d = 2,3,4 or 5 in GENTLE /4fan 
for calculations of Born cross sections, QED FF 
corrected cross sections, QED SF corrected cross 
sections, QED SF corrected angular distributions. 
A calculation time below 15 min seems to be pos- 
sible for d = 3 for all solvable processes, and for 



d = 4 for selected processes only. 

4. The semi-analytical approach 

SA calculations result to nice analytical formu- 
lae, which depend only on the input parameters 
and on the remaining integration variables. These 
formulae show symmetries, which allow a deeper 
physical understanding of the underlying process. 

Specific four fermion processes are given by the 
imaginary part of three-loop diagrams. There- 
fore, some formulae obtained in the SA ap- 
proach agree with results of three-loop calcula- 
tions. However, one has to admit flavour identifi- 
cation in the four fermion final states. Therefore, 
one is not summing over complete weak multi- 
plets. Furthermore, one is not interested in a 
complete analytical phase space integration be- 
cause one wants to calculate distributions and ap- 
ply cuts. 

In contrast to the MC, the SA approach can 
distinguish analytical zeros from contributions, 
which are multiplied by (very) small factors. This 
found an interesting application in the discussion 



of the interferences in SM Higgs production [ 23 



and can be extended to SUSY Higgs production. 

4.1. Cross sections, distributions and cuts 

So far, the SA approach [ g, |J, || || uses the 
following parametrization of the phase space ^| 



dn 



Pl i=l 



(11) 



\A(s,si,s 2 ) yj\{si,m\,ml) 



8.s 



y/Ks 2 ,ml,ml) 



8si 

dsidszdilodili dQ 2 , 



k\ and k 2 are the four-momenta of the initial elec- 
tron and positron. The fermions fi in the final 
state have four-momenta pi and masses m%. The 
invariants s, s\, and s 2 are 



s = (fci+fc 2 ) 2 , 



si 



(P1+P2) 2 , s 2 = {p 3 +p 4 f 



(12) 



3 Other paramctrizations are considered for different 
topologies of Feynman diagrams. 
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One has to substitute cIQq = 2nd cos 9 = 2itdcQ 
to integrate over the rotation angle around the 
beam axis. 9q is the angle between the vectors 
(pi +P2) and k\. The spherical angles of the mo- 
menta pi and p*2 {p*3 and pi) in their rest frames 
are fii (^2): ddi = d cos Oidifii = dcidipi. The 
kinematical ranges of the integration variables 
are: 

(mi +m 2 ) 2 < s± < (y/s — m 3 — m 4 ) 2 , 
(m 3 +TO 4 ) 2 < s 2 < (y/s - % /sT) 2 , 
-1 < d< 1, 
< 0i < 2tt, i= 0,1,2. (13) 

Usually, the squared matrix clement is inte- 
grated over all six angles Ci,(f>i. The remaining 
two integrations are done numerically. As a re- 
sult, the cross section is obtained in the following 
form: 



a(s) 



ds-. 



dsi 



ds\ds2 



(14) 



The boundaries of the integrations s^ 2 allow the 
implementation of cuts in s\ and S2- These 
cuts are very important to separate intermedi- 
ate photons, Z's, W's and Higgs bosons in the 
diagrams lc and If from the background. The 
double differential cross section factorizes, 



d 2 a 
ds\ds2 



i 

xGi(s, si, s 2 ), 



(15) 



where Ci(. . .) are functions depending on the cou- 
plings and the invariants and Gi(. . .) are kine- 
matical functions resulting from the six-fold an- 
alytical integration. They depend on the invari- 
ants only. The summation runs over the interfer- 
ences of different topologies. The double differ- 
ential cross sections allow the calculation of the 
distributions 

da da da , . 

dsi ds2 dhj\2 

where E12 is the sum of the energies of the parti- 
cles /1 and / 2 , 



E\2 = Pw + P20 



si - s 2 



2yfs 



(17) 



Cuts are possible on all parameters in equation 
( |l6| ) simultaneously. 

To be more differential, the integration over Cq 
can be left for numerical integration [ |2(| . The 
results are triple differential cross sections, 



d 3 a 



ds\ds2dcQ 



^2 Ci ( e > A' /2> h, U, s, 81, s 2 ) 

i 

xGi(s,si,S2,c ) (18) 



with the remaining numerical integrations 



a(s) 



asids 2 dco 



The parameters si, S2 and cq allow the construc- 
tion of additional distributions, 



(20) 



da da da 
dc ' dpj 2 ' dy 12 ' 

where p\ 2 is the sum of the transversal momenta 
of the particles /1 and / 2 against the beam axis, 
and y\2 is the rapidity related to p\ 2 - The new 
parameters depend on si, S2 and Co only, 



T 
Pl2 



yjE 2 12 -( mi +m2) 2 ^l-4 



tanhyi2 = cqJI- (mi + m 2 ) 2 /E 2 



(21) 



Again, cuts arc possible on all parameters in 
equations ( |l6|) and (|20| ) simultaneously. 

All distributions listed in equations (^) and 
( ^ii| ) can be calculated with initial state radia- 
tion too. This is possible because Si and S2 are 
invariants, and Co can be reconstructed after a 
boost due to photon radiation from the initial 
state without use of f£i and f2 2 . 

Alternatively, one of the parameters ci or C2 
can be left for numerical integration. The corre- 
sponding cross section formulae are identical to 
( |l8| ) and ( |l9| ) with cq substituted by c\ or 02- No 
new functions C;(. . .) and Gi{. . .) appear in the 
case of NC32 processes. The distributions 



da 



da da 



da 



dpio ' dpf ' dy t ' dci 2 ' 
are now calculable because 



1,2 



(22) 



P10 



(si + m 2 — m 2 ){s + Si 



S2j 



4s ly /s 
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-ci 



v / A(si,TO^,m|)A(s, si,s 2 ) 



4si 



T 
Pi 



tanh yi 

C12 



Pw 



T 



ci 



yjl-ml/p 2 w and (23) 

1 - (1 -cf)A(s,5i,a 2 )/(4ssi) 
l + (l-cf)A(s,si,s 2 )/(4ssi)' 



depend on si, S2 and ci only, p\ is the transver- 
sal momentum against the p\ + p2 axis, ?/i is the 
rapidity related to pf , and c\i is the cosine of the 
angle between p\ and pi- For simplicity, the for- 
mula for C12 is given only in the massless limit. 
As in the case with open cq, cuts are possible to 
the parameters in equations (^) and ( p2~|) simul- 
taneously. Similar formulae can be written for the 
case where c 2 is left for numerical integration. 

The distributions ( p3| ) can be calculated in the 
SA approach only at the Born level. The reason is 
that ci(c2) cannot be reconstructed after a boost 
of the subsystem due to photon radiation without 
knowledge of cq and <^>i(</> 2 ). Unfortunately, these 
angles are already integrated out. 

4.2. Limitations and perspectives 

Presently, only final states marked in bold in 
tables 1 and 2 (NC32, CC11) are treated in the 
SA approach. Final states printed in the tables 
in roman (NC48, CC20) are under investigation 
The NC43, NC19 (CC43, CC19) can cer- 
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tainly be treated. 

All processes with identical particles in the final 
state (NC 4- 36) given in table 2 in typewriter 
cannot be treated with cuts on s\ = (P1+P2) 2 and 
s 2 = (pa +Pi) 2 - These cuts imply the same cuts 
on s\ = (pi + P4) 2 and s 2 = (^3 + P2) 2 because 
the identical final particles are indistinguishable. 
They have to be calculated with a phase space 
parametrization, which has the four invariants 
si, S2,si,S2 as integration parameters. Unfortu- 
nately, the transformation to these parameters is 
described by general polynomials of fourth power, 
which have to be inverted. 

The four jet cross section with two quarks and 
two gluons in the final state [ |2^] has to be added 
incoherently to the NC32 and CC11 processes. 
All soft singularities can be eliminated by a cut 



on the invariants si and s 2 . The collinear sin- 
gularities are regularized by a finite quark mass 
m q . However, the differential cross section is en- 
hanced by a factor s/rri^ for every gluon, which 
is collinear with a quark. Experimentally, one 
can discover a four jet event only, if the jets 
are separated by a minimal angle. Such an an- 
gular cut is not possible in the SA approach. 
This shows another limitation of this calculation 
scheme. The problem arises already in four quark 
final final states but without collinear enhance- 
ment. 

Finally, I would like to add some remarks about 
possible future SA calculations. 

The inclusion of anomalous couplings is possi- 
ble for the CC11 process [ The treatment of 
SUSY Higgs cross section is possible everywhere, 
where the corresponding topology with SM gauge 
bosons can be calculated. This includes the s 
channel Higgs diagrams in muon colliders. Four 
fcrmion final states in e~e~ collisions can be cal- 
culated, if the corresponding topologies are calcu- 
lable in e + e~ collisions. 76(77) collisions demand 
one more integration for every photon in the ini- 
tial state. This complication could be compen- 
sated by the missing QED corrections from the 
initial photon. The inclusion of SUSY processes 
with four fermions in the final states is interest- 
ing. Unfortunately, many SUSY processes lead to 
six fermion final states. 

To summarize this section, we emphasize that 
the SA and MC approach are complementary. 
Both approaches have advantages and disadvan- 
tages. Both are needed at different places in a 
future data analysis. 

5. Summary 

The physics of four fermion final states is much 
richer than that of two fermion final states. The 
calculation of four fermion cross sections with a 
theoretical error below 0.5% is challenging. 

In general, complicated kinematical cuts are re- 
quired already for Born cross sections to ensure 
the applicability of perturbation theory and to 
separate interesting events from the background. 
Radiative corrections are necessary. The codes 
used in fits to future data must be fast and accu- 
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rate. 

Not all these features can be met by a code 
based on one calculation scheme only. Flexible 
Monte Carlo programs and fast semi-analytical 
programs together could fulfill all requirements. 

Present codes probably meet the experimental 
precision expected at LEP2 this year. Theoreti- 
cal errors certainly must be reduced to meet the 
accuracy of the LEP 2 data in the future. 
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